The Ohio State University 


AN E-PLANE ANALYSIS OF APERTURE -MATCHED HORN ANTENNAS 
USING THE MOMENT METHOD AND THE UNIFORM 
GEOMETRICAL THEORY OF DIFFRACTION 


D.J. Heedy 
W.D. Burnside 


The Ohio State University 

ElectroScience Laboratory 

Deportment of Electrical Engineering 
Columbus, Ohio 43212 

(NASA-CR- 176841 ) AN E-ELANE ANALYSIS OF N86-26464 

APERT UBE-MATCEEC HORN ANTENNAS USING THE 
MOMENT METHOD AND THE UNIFORM GEOMETRICAL 

THEORY OF DIFFRACTION (Ohio State Univ.) Onclas 

114 p HC A06/MF A01 CSCL 20N G3/32 43425 

Technical Report 716148-1 
Grant No. NSG 1613 
September 1984 


National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


NOTICES 


When Government drawings, specifications, or other data are 
used for any purpose other than in connection with a definitely 
related Government procurement operation, the United States 
Government thereby incurs no responsibility nor any obligation 
whatsoever, and the fact that the Government may have formulated, 
furnished, or in any way supplied the said drawings, specifications, 
or other data, is not to be regarded by implication or otherwise as 
in any manner licensing the holder or any other person or corporation, 
or conveying any rights or permission to manufacture, use, or sell 
any patented invention that may in any way be related thereto. 


30272 -IQ1 | •_ 

REPORT DOCUMENTATION I report no. 2 - 

page [ 716148-1 

4. Title end Subtitle 

AN E-PLANE ANALYSIS OF APERTURE-MATCHED HORNANTENNAS USING TH 
MOMENT METHOD AND THE UNIFORM GEOMETRICAL THEORY OF DIFFRACTIO 

7. Author(s) 

D.J. Heedy, W.D. Burnside 

9. Performing Organization Name and Address 

The Ohio State University 
ElectroScience Laboratory 
1320 Kinnear Road 
Columbus, Ohio 43212 

12. Sponsoring Organization Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 

15. Supplementary Notes 


16. Abstract (Limit: 200 words) 


The moment method and the uniform geometrical theory of diffraction are utilized to 
obtain two separate solutions for the E-plane field pattern of an aperture-matched 
horn antenna. This particular horn antenna consists of a standard pyramidal horn 
with the folloiwng modifications: a rolled edge section attached to the aperture 

edges and a curved throat section. The resulting geometry provides significantly 
better performance in terms of the pattern, impedance, and frequency characteristics 
than normally obtainable. The moment method is used to calcualte the E-plane pattern 
and BSWR of the antenna. However, at higher frequencies, large amounts of 
computation time are required. The uniform geometrical theory of diffraction 
provides a quick and efficient high frequency solution for the E-plane field pattern. 
In fact, the uniform geometrical theory of diffraction may be used to initially 
design the antenna; then, the moment method may be applied to "fine tune" the design. 
This procedure has been successfully applied to a comapct range feed design. 


17. Document Analysis ft. Descriptor! 


b. I dent I flea/ Open- Ended Terms 


C. COSATI Field/Group 


18. Availability Statement 

19. Security Class (This Report) 

21. No. of Pages 


UNCLASSIFIED 

105 



22. Price 


(See ANSI— Z39.18) See |„.t fuc tion« on Reverse OPTIONAL FORM 272 (4-77) 


(formerly NTIS-3S) 
Department of Commerce 



1 






TABLE OF CONTENTS 


Page 

ACKNOWLEDGMENT ii 

LIST OF TABLES iv 

LIST OF FIGURES v 

Chapter 

I. INTRODUCTION 1 

II. BACKGROUND AND THEORY FOR MOMENT METHOD SOLUTION 7 

III. BACKGROUND FOR UNIFORM GEOMETRICAL THEORY OF 

DIFFRACTION 25 

IV. RESULTS AND ANALYSIS 46 

V. SUMMARY AND CONCLUSIONS 84 

APPENDIX 

A. EVALUATION OF THE DIFFRACTION COEFFICIENT FOR A 

SOURCE MOUNTED ON A PERFECTLY CONDUCTING PLANAR 
SURFACE SMOOTHLY TERMINATED BY A CIRCULAR 
CYLINDER 87 

REFERENCES 104 


iii 


PRECEDING PAGE BLANK NOT FILMED 



LIST OF TABLES 


Table 

1 . 


MEASURED VSWR VERSUS FREQUENCY FOR COMPACT RANGE 
APERTURE-MATCHED HORN ANTENNA WITH ABSORBING 
MATERIAL ON ROLLED EDGES 


Page 


77 


iv 



LIST OF FIGURES 


Figure Page 

1 . Pyramidal horn geometry. 3 

2. Aperture-matched horn geometry (E-plane side view). 5 

■f *► > 

3. The source ( J-,- ,M-j ) generates the field (E,H) with 

scatterer. 8 

4. The interior field vanishes when the currents (J S ,M S ) 

are introduced on the surface of the scatterer. 8 

-V- -V 

5. The exterior scattered field may be generated by (J S ,M S ) 

in free space. 8 

6 . An electric test source J m is positioned in the 

interior of the scattering region. 10 

7 . A planar strip dipole with edges at xj and X 3 and 

terminals at X 2 . 12 

8 . The current-density distribution J on the 

sinusoidal strip dipole. 12 

9. Nonplanar strip dipole with edges at s^ and t^ and 

terminals at 0. 14 

10. An electric strip monopole and the coordinate system. 15 

11. Perfectly conducting polygon cylinder with parallel 

magnetic line source . 17 

12. Electric test probes 1 and 2 are moved to the 

conducting surface. 17 

13. Two-dimensional pyramidal horn geometry for moment method 

solution. 19 

14. Two-dimensional pyramidal horn with curved throat geometry 

for moment method solution. 20 

15. Two-dimensional pyramidal horn with rolled edge geometry 

for moment method solution. 21 


v 



Figure 


Page 


16. Two-dimensional aperture-matched horn geometry for moment 

method solution. 22 

17. Two-dimensional geometry for a source mounted on a 
perfectly conducting planar surface smoothly terminated by 

a circular cylinder. 26 

18. Two-dimensional geometry for pyramidal horn with rolled edge 

for UTD solution. 26 

19. Geometry for a perfectly conducting two-dimensional 

surface with a discontinuity in surface curvature. 28 

20. Two-dimensional geometry for the lit region (0 <<|><tt) 
field for a line source mounted on a perfectly conducting 
planar surface smoothly terminated by a circular cylinder. 31 

21. Two-dimensional geometry for the shadow region ( tt<<j><2tt) 
field for a line source mounted on a perfectly conducting 
planar surface smoothly terminated by a circular cylinder. 32 

22. Two-dimensional aperture-matched horn geometry for UTD 

solution. 36 

23. Circular aperture field distribution. 37 

24. Two-dimensional wedge geometry for calculating the incident 

field from the throat. 40 

25. Two-dimensional aperture-matched horn geometry for UTD 

solution with a modal expansion of the throat field. 44 

26. Moment method solution for standard pyramidal horn. 47 

27. Moment method solution for pyramidal horn with rolled edge. 49 

28. Moment method solution for aperture-matched horn. 52 

29. Moment method and UTD solutions for aperture-matched 

horn as the frequency (f) is varied. 54 

30. Moment method and UTD solutions for aperture-matched horn 

as the radius of curvature (A^) of the throat section is 
varied. 59 


vi 



Figure 


Page 


31. Aperture-matched horn geometry for a large value of the 

radius of curvature (AjJ of the throat section. 63 

32. Moment method and UTD solutions for aperture-matched horn 
as the radius of curvature (A 2 ) of the rolled edge section 

is varied. 64 

33. Moment method and UTD solutions for aperture-matched horn 

as the half-flare angle (a) is varied. 67 

34. Moment method and UTD solutions for aperture-matched horn as 

the aperture width (h) is varied. 72 

35. Compact range aperture-matched horn geometry (E-plane 

side view). 75 

36. Comparison of calculated and measured E-plane field patterns 

for compact range aperture-matched horn. 78 

37. Measured phase versus azimuth angle for compact range 
aperture-matched horn with absorbing material on rolled 

edges. 83 

A.l. Comparison of UTD diffraction coefficient and numerically 
derived diffraction coefficient of Burnside and Chuang as 
the radius of curvature (a 2 ) of the cylinder is varied. 95 

A. 2. Comparison of UTD diffraction coefficient and numerically 
derived diffraction coefficient of Burnside and Chuang as 
the source distance ( p* ) is varied. 100 


vii 



CHAPTER I 


INTRODUCTION 


The use of horn antennas for feeds in compact range measurement 
systems and other applications creates a need for an efficient means of 
designing such antennas. Often it is desirable to do as much of the 
design work as possible "on paper" before constructing a physical model 
to test. Such an approach is usually cheaper and less time-consuming. 

One may make use of one or more analytical techniques to obtain a 
fairly accurate idea of how closely the particular antenna design will 
meet the required specifications before building it. For example, 
aperture integration may be used to find the main beam of the field 
pattern but it requires a knowledge of the field distribution across the 
aperture. Unfortunately, the exact distribution cannot be determined 
without a physical model so an approximate distribution, such as the 
geometrical optics field, is normally used. 

The moment method can provide an "exact" solution to the problem. 

An integral equation which the fields must satisfy is found. From this 
equation, a system of N simultaneous linear equations in N unknowns is 
obtained and solved using linear algebra techniques. Typically, these 
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N unknowns specify the current distribution associated with the body, 
and the result may be used to calculate the fields. For bodies which 
are large with respect to wavelength, the number of unknowns is large, 
and the computation time necessary to obtain a solution may become 
intolerable. Thus, the moment method is an excellent low frequency 
method. It can be used to analyze electrically large structures; 
however, it may require large amounts of computation time at the higher 
frequencies. 

The uniform geometrical theory of diffraction is a good high 
frequency method which has the added advantage that it breaks down 
"gracefully" as the frequency decreases. This method modifies the 
geometrical optics solution by including a diffracted field such that 
the total field is smooth and continuous. A diffracted field may be 
associated with a discontinuity in the incident or reflected field as 
well as with a discontinuity in the rate of change or "slope" of the 
fields [1,2,3], These discontinuities in the fields are the result of 
discontinuities in the extent or curvature of perfectly conducting 
bodies. 

The uniform geometrical theory of diffraction (herewith referred to 
as UTD) can in many cases provide accurate results with relatively short 
computation time. For instance, it has been used to calculate both 
E-plane and H-plane patterns for the standard pyramidal horn [4,5,6] 
such as shown in Figure 1. 
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Figure 1. Pyramidal horn geometry. 
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In this report, an aperture-matched horn antenna with the geometry 
shown in Figure 2 is considered. In comparison to the standard 
pyramidal horn, the aperture-matched horn has a smoothly curved surface 
attached to the aperture edges. In this case the curved section is 
circular in shape, i.e., the radius of curvature is constant. Another 
curved section is added between the flared walls of the horn and 
waveguide junction. Thus, the sharp edges of the pyramidal horn are 
eliminated. The rolled edges allow much of the energy striking the edge 
to travel along the curved surface gradually shedding away from the 
surface into the sidelobe and backlobe regions. The result is a broader 
main beam with lower sidelobes and backlobes. In addition, the curved 
surfaces decrease the amount of energy reflected back into the throat of 
the horn resulting in a lower VSWR (voltage standing wave ratio). This 
property is important for feed horns in compact range systems where the 
antenna is used to transmit and/or receive. 

The method of moments and the UTD are used to calculate the far 
field E-plane pattern of the aperture-matched horn antenna. The moment 
method may be conveniently used to find both E-plane and H-plane 
patterns at low frequencies when the horn dimensions are small with 
respect to wavelength. It may also be used at higher frequencies, but 
the computation time increases. The UTD solution discussed in this 
report provides good results for the E-plane pattern at higher 
frequencies with much less computation time. Thus, the pattern 
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behavior can be examined over a specified bandwidth. Also, the moment 
method solution is used to calculate the VSWR of the antenna. 

Chapter II discusses the moment method solution for the E-plane 
pattern while Chapter III discusses the UTD solution for the E-plane 
pattern. In Chapter IV, a number of geometries are considered. The 
moment method and UTD solutions for the E-plane are compared and a brief 
analysis of the effects of changing various parameters of the horn 
geometry on the pattern and VSWR is made. Chapter V is a short summary 
with appropriate conclusions. 
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CHAPTER II 


BACKGROUND AND THEORY FOR MOMENT METHOD SOLUTION 


A moment method solution was obtained by using the reaction concept 
of Rumsey [7], Consider the exterior scattering problem illustrated in 

4 4 

Figures 3, 4, and 5. Impressed electric and magnetic currents (Ji »M^ ) 

4 4 “ 

generate electric and magnetic field intensities (E,H) in the presence 
of a conducting body in free space. 

From the surface-equivalence theorem of Schellkunoff [8], an 
equivalent problem may be obtained by replacing the body by the 
following surface current densities: 

■* A. > 


and 

M s = E x n (2.2) 

A 

where the unit vector n is the outward directed normal to the surface S. 

4 * *4 

By definition, the source currents (J.. ,M. ) generate the incident fields 

4 4 

(E..,H^) in free space. The scattered fields are defined as 
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Figure 3. 


-* * ■* 

The source (Ji ,M-j ) generates the field (E,H) with 
scatterer. 



■> 

Figure 4. The interior field vanishes when the currents (J s » 
are introduced on the surface of the scatterer. 


I \ 

I f-E’.-H* ) \ 

\ \ 
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Figure 5. The exterior scattered field may be generated by (J S ,M S ) 
in free space. 


2 + 



(2.3) 


* 




and 



(2.4) 


* + 

The surface currents generate the scattered fields (E S ,H S ) exterior to 
the body and (-E-j,-H.j) in the interior region. 

One now places an electric test source J m in this region as 
illustrated in Figure 6. Since there is a null field in the interior, 
the reaction of this test source with the fields of the other sources is 
zero. By reciprocity, this reaction is equivalent to the reaction of 
the remaining sources with the field of the test source such that 


//( J •£ - M »H )ds + - M. *H )dv - 0 

1 1 v s m s m i m l m 


(2.5) 


■* -► 

where (E m ,H m ) are the free space fields of the test source. This 
equation is a statement of the "zero-reaction theorem" of Rumsey [7] and 
is used to solve this problem. 

Let us determine the surface current distributions (J S ,M S ) from 

which the scattered fields may be calculated. To do this, one expands 

these functions in finite series with N unknown expansion coefficients. 

For our purposes, it is assumed that the body is a perfect conductor of 

electricity so that M s vanishes. Let us also consider a two-dimensional 
■> 

problem in which J s is only a function of the position l around the 


9 




Figure 6. An electric test source J m is positioned in the interior of 
the scattering region. 


contour C of a cylinder. For the TE polarization in which the electric 

* 

field is transverse to the z-axis, a magnetic line source is used and J-j 
is zero. The integral equation. Equation (2.5), becomes 



•E dl = // M.«H ds 
m i m 


( 2 . 6 ) 


Let us represent the electric current distribution as 


o s u) = ! v„(*) 


(2.7) 


where the complex constants (I n ) are samples of the function J s (t), and 

■> 

the vector basis functions are denoted by J n (t). The test source and 
the basis functions have unit current density at their terminals. 
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Substituting Equation (2.7) into Equation (2.6), one obtains the 
following simultaneous linear equations: 


N 

l I Z = V with m = 1,2,3, ...N 
n mn m 


where 


Z 

mn 


= -/ J (l) 
J n n 


E dt 
m 


= -iju) 

m m 



(2.9) 


V = 
m 



Hds = J J U)* E.di 
m mm l 


( 2 . 10 ) 


and the integrations extend over the region where the integrand is 
non-zero. 

For computational speed and storage, it is advantageous to have a 

symmetric impedance matrix Z mn . Furthermore, the test sources should 

yield a well -conditioned set of simultaneous linear equations. For 

these reasons, and to obtain closed forms for some of the integrals in 

Equations (2.9) and (2.10), one uses test sources J m of the same size, 

* 

shape, and functional form as the expansion functions J n . Finally, the 
interior test sources are placed a small distance 6 from the surface S. 
and the limiting form of the integrals is taken as <S tends to zero. 

For ease of computation, the basis functions are chosen to be 
sinusoidal strip dipoles. Such a planar strip dipole is illustrated in 
Figures 7 and 8. This dipole lies in the xz-plane and has infinite 
length in the z-direction. The surface-current density is 
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y 




Figure 7 


y 


L 


Figure 8 




J 


*2 



X 


A planar strip dipole with edges at xi and X3 and 
terminals at X£. 



. The current-density distribution J on the sinusoidal strip 
dipole. 
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( 2 . 11 ) 


♦ a sin(k(x-x 1 )) 

J ” X sin(k(x2-x 1 )) 

for x 1 < x < x^, and 

> a sln(k(x 3 -x)) 

J ” x sin(k(x 3 -X2)l 

for X 2 < x < x 3> Figure 9 illustrates a strip V-dipole. The distance 
along the dipole arms is measured by the coordinates s and t, 
respectively, with the origin at the terminals 0. The surface-current 
density is 

» . s1n(k(sj-s)) 

J = ' s (2-13) 

on arm s, and 

♦ . sin(k(tj-t)) 

J=t slnfktj) < 2 - 14 ) 

A A 

on arm t. The unit vectors s and t are perpendicular to the z-axis. In 
both dipoles, the current density vanishes at the end points and is 
unity at the terminals 0. The current density is continuous across the 
terminals, but it does have a slope discontinuity there. Note that when 
= 180°, the V-dipole reduces to the planar dipole. 
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Figure 9. Nonplanar strip dipole with edges at si and ti and 
terminals at 0. 
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The fields of the sinusoidal strip dipole may be obtained by 
considering the superposition of two strip monopoles with sinusoidal 
current distributions and a common end point. The fields of the strip 
monopole shown in Figure 10 are discussed in reference [9], Using these 
results, let us calculate the various elements Z mn of the impedance 
matrix and the elements V m of the excitation column in Equation (2.8). 

Consider a perfectly conducting polygon cylinder with contour C 
which may be open or closed. Let J S U) denote the surface-current 
density induced on the cylinder. In this case, the cylinder is open, so 
currents will flow on both sides of the thin conducting surface, and J s 
will denote the total current density. 


A 



Figure 10. An electric strip monopole and the coordinate system. 
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Figure 11 illustrates a perfectly conducting polygon cylinder 

■* 

illuminated by a parallel magnetic line source M^ . Let Ij and I 2 denote 
the current density J s at the corners of the polygon. The current 
density Os vanishes at the edges 0 and 3. Let us define two strip 
dipole mode currents on the cylinder. Mode 1 extends from point 0 to 
point 2 and has terminals at point 1. Mode 2 extends from 1 to 3 with 
terminals at 2. Each mode has a sinusoidal current distribution and 
unit terminal current as in Equations (2.13) and (2.14). Now let us 
represent J s (t) as the superposition of the two modal currents with 
weightings Ij and I 2 . This gives a piecewise-sinusoidal expansion for 
J s (fc) with two unknown constants 1^ and I 2 . 

In the exact solution, the tangential electric field vanishes 
everywhere on contour C. Thus, if one moves an electric test probe to 
the conducting surface, as in Figure 12, the open-circuit voltage at 
its terminals will read zero. To determine N current samples, one makes 
N independent probing tests. The probes may be real (thin-wire 
V-di poles) or hypothetical (electric line sources or strip dipoles). 

Now suppose one adjusts the currents I n until all the probes read zero. 
This procedure yields a stationary solution for the currents I n and, 
under favorable conditions, tends to the rigorous solution as N 
increases. 

Let Z mn denote the mutual impedance between test-probe m in Figure 
12 and mode current I n in Figure 11. The open-circuit voltage 
induced in the probe is the sum of the voltage contributions from 
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Figure 11. Perfectly conducting polygon cylinder with parallel 
magnetic line source ft.. 



Figure 12. Electric test probes 1 and 2 are moved to the conducting 
surface. 
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Js and Mi. This voltage must vanish at each probe, leading again to 
Equations (2.8) to (2.10). 

Linear algebra techniques are used to solve these simultaneous 
linear equations and determine the current distribution J s ( a ) and the 
scattered fields (E S ,H S ). The theory and method outlined above, as well 
as the appropriate computer programs, were generously provided by J.H. 
Richmond [9], Duality may be used to obtain a solution for the case of 
TM polarization. 

Figures 13 to 16 illustrate the various horn geometries used for 
the moment method solution. The horn is modelled by a two-dimensional 
open cylinder composed of a number of strips of infinite length and 
finite width. For accurate results, the strips should have widths of no 
more than a quarter of a wavelength. In general, the widths should be 
small enough to obtain a good convergent solution. For the E-plane, the 
source is a magnetic line source located on the end wall of the 
waveguide at the center of the waveguide section. 

Figure 13 is a cross-sectional view in the E-plane of a 
conventional pyramidal horn. The geometry is described by four 
parameters: input waveguide width w, input waveguide length X[_, 

half-flare angle a, and aperture width h. The length X|_ of the input 
waveguide should be at least one wavelength to insure that the 
evanescent modes from the line source have decayed sufficiently. Thus, 
when the electromagnetic waves reach the throat of the horn, only the 
dominant waveguide mode is present. The range of values for the 
half-flare angle a is 0°<a<90°. The walls of the cylinder are 
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Figure 13. Two-dimensional pyramidal horn geometry for moment method 
solution. 

divided into segments which have widths of approximately x/8 on the 
waveguide section and X/5 on the straight flare sections. The segment 
widths are smaller on the waveguide section because they are located 
closer to the line source than those on the flare sections. Thus, the 
fields and currents are stronger on the waveguide section and the 
smaller segments provide a more accurate model of the actual current 
distribution. 

Figure 14 shows a pyramidal horn with a circularly curved throat 
section. Note that the junctions between this curved throat section and 
the waveguide and flare sections are smooth; that is, the waveguide and 
flare walls are attached tangentially to the circularly curved throat 


19 




Figure 14. Two-dimensional pyramidal horn with curved throat 
geometry for moment method solution. 

section. This geometry is described by the same parameters as the 
dihedral horn in addition to the radius of curvature of the throat 
section Aj_. The segments of the throat section have widths of 
approximately x/8 for the same reasons discussed earlier. 

Figure 15 is a pyramidal horn with a rolled edge flare section. 

Note that the rolled edge also has a circular curvature and forms a 180° 
circular arc. The parameter A 2 denotes the radius of curvature of the 
rolled edge section. This section has segment widths of approximately 
A/5. 

Finally, Figure 16 illustrates the most general geometry treated in 
this report: an aperture-matched horn with circularly curved throat and 

rolled edge sections. All six parameters (w, xl, a , h, Aj, and A 2 ) are 

needed to describe this geometry. 
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Figure 16. Two-dimensional aperture-matched horn geometry for moment 
method solution. 
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One advantage of an aperture-matched horn is a relatively low VSWR. 
Only a fraction of the incident wave travelling down the waveguide is 
reflected back to the source. This characteristic makes it desirable 
to calculate and predict the VSWR for such a horn geometry. 

The field in the waveguide section is the sum of an incident wave 
travelling from the source to the aperture and a reflected field 
travelling from the aperture to the source. At any given point in the 
waveguide, the two waves add or subtract depending on their relative 
phase. The result is a standing wave pattern from which the VSWR and 
equivalently the reflection coefficient can be obtained. 

The moment method solution is used to calculate the magnitude of 
the field at a number of points along the waveguide section. The minima 
and maxima of the standing wave pattern are determined from these 
values. The VSWR is given by 


VSWR 


E MAX 

I e min 


MAX' 



(2.15) 


where (E.,* v>H. lllv/ ) and are the maximum and minimum field 

MAX MAX MIN MIN 

magnitudes, respectively. The reflection coefficient p can then be 
calculated from elementary considerations: 


p = 



(2.16) 
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VSWR-1 
|p| = VSWR+1 


(2.17) 


and 


5 = 2kx n - (20-1)* (2.18) 

where k is the wave number and x n is the distance from the nth minima in 
the waveguide to the aperture. 
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CHAPTER III 


BACKGROUND FOR UNIFORM GEOMETRICAL THEORY 
OF DIFFRACTION 

Complete uniform geometrical theory of diffraction solutions for 

the standard pyramidal horn have been obtained for both the E-plane and 

H-plane patterns by Peters, Rudduck, and Yu [5,6], The UTD wedge 

diffraction solutions were utilized in these papers. The corresponding 
two-dimensional geometry for the moment method solution is illustrated 

in Figure 13. 

Burnside and Chuang [10,11] used a hybrid approach to obtain a 
solution for an aperture-matched horn design consisting of a dihedral 
horn with a rolled edge flare section as shown in Figure 15. This 
hybrid approach combined moment method and UTD techniques to obtain a 
numerically derived solution for the diffraction coefficient of a source 
mounted on a perfectly conducting planar surface smoothly terminated by 
a circular cylinder as shown in Figure 17. 

Two such surfaces may be attached as shown in Figure 18 to model 
a horn geometry. A magnetic line source at the vertex generates the 
field from the throat of the horn. Note that the throat region appears 
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meld 

POINT 0 



Figure 17. Two-dimensional geometry for a source mounted on a perfectly 
conducting planar surface smoothly terminated by a circular 
cylinder. 



Figure 18. Two-dimensional geometry for pyramidal horn with rolled edge 
for UTD solution. 
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as an electrically small radiator; as expected, its pattern is smooth 
across the horn flare angle and zero otherwise. Thus, the pattern is 
dominated by three terms: direct throat radiation and two diffraction 

terms from the edge junctions. The resulting solution agrees quite well 
with measurements and complete moment method solutions. 

The aperture-matched horn design shown in Figure 16 has, in 
addition to the rolled edge flare section, a circularly curved throat 
section. Thus, a more accurate model is needed to include the effects 
of the curved section on the source field from the throat. Secondly, to 
decrease computation time a "pure" UTD solution is preferable to the 
hybrid approach. In other words, one needs a UTD solution of the 
diffraction coefficient for a source mounted on a perfectly conducting 
planar surface smoothly terminated by a circular cylinder. 

Such a diffraction coefficient is obtained by empirically modifying 
a previous solution of a diffraction coefficient for a perfectly 
conducting two-dimensional surface with a discontinuity in surface 
curvature [12] as shown in Figure 19. The diffraction coefficient and 
the necessary parameters are given by 

e -j*/4 C h (C 1 )F(X 1 ) - C h U 2 )F(X 2 ) 

D h = cos* + cos $ 1 (3.1) 


where 


C h (5 l,2 ) 


2cos(HP“) 

a l»2 


2 J ttX F(X) + m 1> 2 x/F q (Si >2 ) 


-j*/4 -jU, ? r/12 
e e ’ 


(3.2) 
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REFLECTED RAY 



Figure 19. Geometry for a perfectly conducting two-dimensional surface 
with a discontinuity in surface curvature. 
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(3.3) 


h,2 " 2m l,2 C0s 


,±*Ly 

l 2 ) 


m l ,2 (l/2ka^ 2 ) 1/3 


(3.4) 


L = 


ss' 
s+s 1 


(3.5) 


X = 2kL cos 


( 2 ) 


(3.6) 


1,2 


2 

ka^ 2 (cos<j> + cos <(> ' ) 


4cos (l^L) [1+(I + L.) -lli cos(ilL)l 
2 s s 2 2 


(3.7) 


OO # 2 

F(x) = 2j/x e"* x / e -J< ' dt 

/x 


(3.8) 


and 


q (5 1>2 ) 


1 » v'(x) -jc. 9 r 

*1 ' dT 


(3.9) 


Note that V(t) and W 2 (t) are Fock type Airy functions given by 


2jV(t) = W 1 (t) - W 2 (x) 


(3.10) 
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where 


1 ®-je 3 /3 

“a = * e dt (3 ’ u) 


with e being a positive number however small; further, 

9V 2 (t) 

V(r)= — 

and 


W 2 (t) 


3W 2 (t) 

3t 


(3.12) 


(3.13) 


This diffraction coefficient is derived for the case where the source is 
off the conducting cylinder (<}>'*0). 

Let us assume that a^ becomes infinite and <(>'=0 so that grazing 
incidence occurs as shown in Figures 20 and 21. Notice that there now 
exists a lit region (0<4><ir) and a dark region (n<<t><2ir) where the 
incident field is blocked by the conducting cylinder. The necessary 
modifications are described in Appendix A. The resulting diffracted 
field for TE polarization is given in the lit region (0<<|><ir) by 

g-jks 

H(Q) = zH^(Qj)D^ (3.14) 


where H^Q,) is the field incident on the junction at the point Q., 
J J 

s is the distance from Q to the field point, and 

J 
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FIELD POINT 
Q 


Figure 21. Two-dimensional geometry for the shadow region ( ir<4»<2ir ) 
field for a line source mounted on a perfectly conducting 
planar surface smoothly terminated by a circular cylinder. 




(3.15) 



cos(ilil) F(X X ) + c h u 2 ) F(X 2 ) 

COS(|) + COS^' 


W 



■) 


\ JlK F(X) + m 2 ^ £ q* 

-jir/4 -j(*- 2 ) 3 /12 
e e 


Z 2 = - 2m 2 cos (-g— ) 
m 2 = (1/2 ka 2 ) 1/3 


ss* 

L = s + s' 


9 <b-d> ' 

X = 2kL cos^( 2 ) 


X, = 


2 

kL(cos<(i + cos4>‘ ) 

n ? j, 4>“ 'fr * >. 
8cos (“2 ) 


and, 


Xo = 


2 

ka 2 (cos<|> + cos<}>' ) 

.d>— d) 1 . r ,1 1 v , (b— <b * . .. 

4cos + ( 7 + T" ) a 2 / 2 cos ("V - )! 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 
(3.22) 
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In the shadow region ( Tr<<j><2Tr) , the diffracted field is given by 


H(Q) = zH ( Qj ) C x 


e -jks 


d 


where 


m 2 e 


•jkt 


__ F(X d ) *, % 
2S/T ' q (C) 


- j tt/4 


m 2 = (l/2ka 2 ) 1/3 


t - &2 1 4*2 “ ^1 1 


m 2 t 



d _ k s‘s d K 2 
~ s'+s d 2m| 


h l (q) 

c = H s (Q) 


e -jks _ 



_ e - 
H 1 (Q) + H^Qj) D h ~7s 



e -jks d 
H (Q J )x h ~W 


SB 




(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


(3.29) 
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and is the distance from the creeping wave launch point Q-p to the 
field point. The functions F(x) and q*(x) are defined in Equations 
(3.8) to (3.13). This coefficient is found to give good results for 
s' > 1.5X and > 0.5\ (see Appendix A for details). 

The aperture-matched horn is modelled by the two-dimensional 
geometry shown in Figure 22 with the x-axis being the axis of symmetry 
of the horn. A plane wave with the E field polarized in the y-direction 
is incident in the waveguide. The field strikes the junctions at Qjh 
and Qji 2 resulting in the diffracted fields Hji and Hj 2 * 


Consider a circular aperture of radius (R^) centered at the apex 
(P) of the horn. The radius (R^) is such that the aperture is at the 
junctions (Qj$j and Qj^) between the curved throat sections and the 
walls of the horn. The creeping wave fields which travel around the 
curved section and strike the junctions (Q C1 and Q lco ) will also 
produce diffracted fields. However, it is assumed that the creeping 
waves have attenuated substantially by the time they reach the second 
junction. Thus, to a first-order approximation, this term may be 
neglected. Considering only first-order terms, the field distribution 
across the aperture has three contributions: the incident waveguide 

field H^(Q) and the two diffracted fields Hj^(Q) and H^Q) from 
junctions Qj^ and Qj^» respectively. 

The horn aperture with the various regions are illustrated in 
Figure 23. The aperture fields are given by 


H(Q) = H i (Q) + Hj^Q) + Hj 2 (Q) 


(3.30) 
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in region II, and 


H(Q) > hJj(Q) + H^ 2 (Q) (3.3!) 

in regions I and III. The component fields at the point Q = in 

the aperture are 


H 1 (Q) = z e A 


(3.32) 


^(Q) - 


-jksi 


z H^Q,..) D 


Jlr hi /sj" 


-jksj 


: u i 


in regions II and III 


2 H (Q JU> C 1 T hl 1" region I 


(3.33) 


and 


hJ 2 (Q) ■ 


■jks 2 


z H 1 (Q 110 ) D 


J12' h2 /sj 


- jkS 2 


* i 

_ ui 


in regions I and II 


2 H (Q J12> C 2 T h2 — 7ZT re9i0 " 111 


(3.34) 


where Sj 2> 2 , Cj 2> and 2 are as defined in Equations 

(3.14) to (3.29) applied to junction Q . For plane wave incidence 

J 1 1 yC. 

in the near field, the constant ^ 1S Qiven by 
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(3.35) 


m/2 q (0) e"^ 4 + /ksi 2 (1-2H 1 (Qjh 2 ) ^ 

C = ’ : 1 

X » 2 2m/? q*(0) e-^ 4 - /lcsTT 

where 


m = (l/2kA 1 ) 1/3 


(3.36) 


and q*(0) is the appropriate Fock integral evaluated at zero. 

Once the aperture field distribution is known, one may replace the 
original geometry with the wedge geometry shown in Figure 24. 

Neglecting, for the moment, the finite termination of the perfectly 
conducting wedge at the rolled edges of the horn, the resulting geometry 
is an infinite two-dimensional wedge of angle equal to 2a. 

The magnetic field (Hp must satisfy the scalar Helmholtz 
equation in this source-free wedge-shaped region which is given by . 


2 i 2 i 
v + K = 0 


(3.37) 


where k is the free space wave number. The accompanying boundary 
conditions are obtained from the condition on the tangential electric 
field at the perfectly conducting boundaries such that 


Ep(p»<J>) = 


1 3Hj(p,«(>) 

jo >e = 0 


(3.38) 
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at $ = ± a. Also, ( p , <fr) must satisfy the radiation conditions as 
p • for e Jwt time dependence. Thus, an appropriate solution for 

h’(p,*) is 

00 

H*(p,*) = l A cos ) H< 2) (kp) (3.39) 

where A n are constants to be determined. For kp large, as in the far 
field, one finds that 

H i 2) ' k| ’> : idl J " e ' JkP • (3.40) 

For n=0, this approximation is good to within about 10% error for 
p > 0.2X. Then, Equation (3.37) becomes 


g-jkp 00 

H ^( P* 4>) = /— l K n cos(— <(>) (3.41) 

• n=0 

where 

S n ■ J " A n /IT • (3.42) 

On the circular arc aperture, p * so that 
"JkR. 

0 “ 00 p ^ 

H z (R A >+) = 7 r^ \ cos ( « 4) * H zA U) (3.43) 
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where H zA (<f>) is the known aperture distribution. Equation (3.43) is 

simply a Fourier cosine series representation of H zA (<|>). The 

coefficients K are found by multiplying both sides of Equation (3.43) 
mr 

by cos (~ <|>) and integrating over the interval -a < <|> < a. The 

result is 

e n /R/\ -jkR A a nir 

“ 2a e / H z y\( ♦) cos t~a ♦) d 4> (3.44) 

-a 

where 

1, n = 0 
G n 

2, n * 0 (3.45) 

which is Neumann's number. 

Once the coefficients (A p ) have been calculated. Equation (3.41) 
gives the field (H z ) for large p. Provided H zA (<f>) is a relatively 
smooth function across the circular aperture, the series will converge 
rapidly and only the first few terms or modes will be significant. The 
infinite series of Equation (3.41) becomes a finite series for the 
incident throat field of the horn such that 


Hj(p,*) 



? r ( nir 

I A cos(— 
n=0 n 


40 


(3.46) 


where N is a finite number such that A n for n > N are relatively small. 
Note that each mode in the series has the same p dependence as a line 
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source field coming from the apex, and each can be treated as such in 
the subsequent UTD analysis. 

The geometry of the problem now reduces to that of Figure 25. Note 
that it is identical to the model used by Burnside and Chuang [11] 
except that the incident field from the throat of the horn is modified 
to account for the curved throat sections. In a first-order solution 
there are three contributions: the incident throat field and the two 

diffracted fields from the junctions Qj2i anc * Qj22* The incident throat 
field is given by 


,mr 


H i (Q) = 


\ e-Jkp N 
z > — > ft cos (~ <f>) in region I 

/p n=0 n 


in regions II, III, IV, and V 


(3.47) 


where the A n are defined in Equation (3.44). 
+d2 

and H from the junctions at and Qj^ 


The diffracted fields 
are, respectively, given by 


-jkpi 


z H 1 (Q 101 )D 


J2r hi 


H dl (Q> ■ 


d 

-jkpi 


2 H "WVhi - 7 % 


in regions I and V 

in regions II and III 
in region IV 

(3.48) 
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and 


~J kp 2 


A j ^ \ ^ 

2 H ^J22^h2 Tpf in regions I and II 


5 d2 (Q) - 


d 

-jkp 2 


A -j 

z H'(Q J22 )C 2 T h2 7$~~ in regions IV and V 


in region III 


(3.49) 


where p 1 ? , 2> D hl 2 » ? , and T hl>2 are as defined in Equations 

(3.12) to (3.27) applied to the Q J21 2 junction. For the far field for 
a line source incident on the junction one finds that 


m/2 - q*(0) e" jll/4 - SW; 

Cl » 2 " 2m rr q *(o) e-J ^ 4 - m; 


(3.50) 


where 


m = (1/2 k A 2 ) 1/3 


(3.51) 


and q*(0) is the appropriate Fock integral evaluated at zero. 
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CHAPTER IV 


RESULTS AND ANALYSIS 

The purpose of this chapter is twofold. Results obtained from the 
theory discussed in the previous two chapters are presented and analyzed 
to determine the validity of the theory. Secondly, the value of the 
particular aperture-matched horn design considered in this report as 
well as the effect of changing various design parameters are discussed. 

Consider first the two-dimensional geometry for the standard 
pyramidal horn shown in Figure 26(a). The frequency is 10 Ghz, and the 
remaining parameters are given in the figure. The pattern obtained by 
the moment method solution is shown in Figure 26(b), and the calculated 
VSWR is 1.35. 

Next, a rolled edge of radius = 7.62 cm is added to this 
pyramidal horn as illustrated in Figure 27(a). Note that this is the 
two-dimensional geometry treated by Burnside and Chuang in Reference 
[11], The resulting moment method pattern is shown in Figure 27(b). As 
expected, the main beam of the pattern is broader than that of the 
pyramidal horn, and the side lobes are much lower. The sharp edges of 
the pyramidal horn diffract the energy into the main beam region where 
it interferes with the incident field from the throat. The rolled 
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(a). Two-dimensional geometry for pyramidal horn (E-plane), 
f = 10 GHz, w = 1.016 cm, X|_ = 10.16 cm, a = 45°, 
h = 7.85 cm. 

Figure 26. Moment method solution for standard pyramidal horn. 


47 



o 

o 



(b). E-plane field pattern for pyramidal horn obtained by 
moment method solution. VSWR = 1.35. 


Figure 26. (Continued). 
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(a). Two dimensional geometry for pyramidal horn with rolled edge 
section (E-plane). 

f = 10 GHz, w = 1.016 cm, x L = 10.16 cm, a = 45°, 
h = 7.85 cm, A 2 = 7.62 cm. 


Figure 27. Moment method solution for pyramidal horn with rolled edge. 
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(b) E-plane field pattern for pyramidal horn with rolled edge 
section obtained by moment method solution. VSWR = 1.35. 

Figure 27. (Continued). 
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edges, on the other hand, allow the energy to creep around the curved 
section and gradually shed off into the side. lobe region. The 
calculated VSWR is 1.35. 

Figure 28(a) illustrates the aperture-matched horn with a curved 
throat section with a radius of 3.05 cm. The corresponding pattern 
shown in Figure 28(b) is almost identical to that of Figure 27(b). 
However, the calculated VSWR is 1.04. This is almost an 87% reduction 
in the reflection coefficient at the throat of the horn. Thus, the 
rolled edges improve the pattern by making it broader and reducing the 
side lobes; in addition, the curved throat section serves mainly to 
reduce the VSWR. 

The frequency dependence of the field pattern and VSWR of the 
aperture-matched horn is examined in Figure 29 using the moment method 
approach. At the low frequencies, the antenna is small with respect to 
a wavelength. As is typical for electrically small antennas, the 
pattern has a broad main beam with relatively high side lobes. As the 
frequency increases, the main beam becomes narrower, and the side lobes 
become lower. The VSWR decreases as the horn, particularly the radius 
of curvature of the throat, becomes larger with respect to the 
wavelength. 

The results of the UTD approach are also shown in Figure 29. UTD 
is a high frequency technique, and the field plots clearly show this 
property. The solution, however, agrees quite well with the more exact 
moment method results at frequencies as low as 6 GHz (X = 5 cm). At 
this frequency, the radius Aj of the throat section is about 0.6X. This 
is quite close to the low frequency limit of the diffraction coefficient 
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(a). Two-dimensional geometry for aperture-matched horn 
(E-plane). 

f = 10 GHz, w = 1.016 cm, X 2 = 10.16 cm, a = 45°, 
h = 7.85 cm, = 3.05 cm, A 2 = 7.62 cm. 


Figure 28. Moment method solution for aperture-matched horn 
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(b). E-plane field pattern for aperture-matched horn obtained 
by moment method solution. VSWR = 1.04. 


Figure 28. (Continued) 




(a). Two-dimensional geometry for aperture-matched horn (E-plane), 
w = 1.016 cm, a = 45°, h = 7.85 cm, = 3.05 cm, 

- 7.62 cm, f varies. 



(b). E-plane field patterns for f = 2 Ghz. VSWR = 2.52. 


Figure 29. Moment method and UTD solutions for aperture-matched horn 
as the frequency (f) is varied. 
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•(g). 


(h). 

Figure 29. 
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E-plane field patterns for f = 14 Ghz 
(Continued) . 
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VSWR = 1.02. 
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Dpc which is valid for Aj > 0.5X. In addition, the apparent source 
distance s' (denoted by R s in Figure 25) for the rolled edge is only 
1.1X. This is outside the range of validity of DpQ, and it probably 
accounts for the disagreement with the moment method solution in the 
region outside the flare angle. Above 6 GHz, the UTD patterns are 
within 1 dB of the moment method patterns in the main beam region. 

There are several geometrical parameters of the aperture-matched 
horn which may affect the performance of a particular design. For 
instance. Figure 30 illustrates the pattern behavior as the radius of 
curvature Aj of the throat section is varied. The results are obtained 
from the moment method. As Aj decreases, the geometry approaches that 
of Figure 27(a). In fact, the pattern in Figure 30 (A} = 0.75 cm = 

0. 25X) is almost identical to the pattern in Figure 27(b). As A^ 
increases, the VSWR decreases. This is expected since a larger radius 
of curvature provides a more gradual transition from the waveguide to 
horn modes. Thus, there is a smaller reflection coefficient associated 
with the junction. 

The corresponding patterns obtained from the UTD approach are also 
shown in Figure 30. For A^ = 0.75 cm = 0.25X, there is a slight 
disagreement with the moment method pattern which dips at 0°. This is 
probably due to the fact that A^ is only a quarter of a wavelength, 

1. e., the low frequency case. However, another discrepancy occurs for 
larger values of A]_. Specifically, for Aj = 6 cm = 2X, the main beam 
begins to become too broad in the UTD solution. 
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(b). E-plane field patterns for A] = 0.75 cm. VSWR = 1.23. 


Figure 30. Moment method and UTD solutions for aperture-matched horn 
as the radius of curvature (Aj) of the throat section is 
varied. 


59 



ANCLE IN DEGREES 


(c). E-plane field patterns for A^ = 1.50 cm, 



(d). E-plane field patterns for Aj = 3.05 cm 


Figure 30. (Continued). 
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(e). E-plane field patterns for Ai = 6.00 cm. VSWR = 1.03. 



(f). E-plane field patterns for Ai = 9.00 cm. VSWR = 1.02. 


Figure 30. (Continued) 



One possible explanation is suggested by the geometry shown in 
Figure 31. As Aj increases with the other parameters remaining fixed, 
the length R|_ of the straight sides of the horn decreases. The form of 
the modal expansion of the field across the circularly curved throat 
aperture assumes the region is bounded by an infinite two-dimensional 
wedge of angle 2a. However, as Rl becomes smaller, this approximation 
becomes less valid as the wedge structure vanishes. Also, note that as 
Ai increases, the radius of curvature R^ of the throat aperture 
increases. Thus, the throat aperture is closer to the horn aperture. 

The diffracted fields from the rolled edge junctions may then need to be 
included in the field distribution across the throat aperture. 

The effect of varying the radius of curvature A 2 of the rolled edge 
section is examined in Figure 32 using the moment method solution. The 
most noticeable change in the pattern as A 2 increases is a decrease in 
the level of the side lobes. In addition, the VSWR drops slightly. The 
UTD solution results are also shown in Figure 32 and agree well with the 
moment method solution for A 2 greater than one wavelength. However for 
smaller values of A 2 , the UTD patterns exhibit a slight slope 
discontinuity at the shadow boundaries at 45° and 315°. The main beam 
of the UTD pattern disagrees with the moment method results, too. This 
behavior may be attributed to the low frequency limitations of the UTD 
diffraction coefficients. 

Figure 33 illustrates the dependence of the pattern on the 
half-flare angle using the moment method approach. Clearly, as a 
increases, the width of the main beam increases. The calculated VSWR 
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f = 10 GHz 



w = 1.02cm 
a = 45° 

h * 7.85cm 
A,* 6 cm 
A 2 * 7.62 cm 
R l *0.78X 

R a = I.07X 


Figure 31. Aperture-matched horn geometry for a large value of the 
radius of curvature (Aj) of the throat section. 
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(a). Two-dimensional geometry for aperture-matched horn (E-plane), 
f = 10 GHz, w = 1.016 cm, o = 45°, h = 7.85 cm. A} * 3.05 cm, 
A 2 varies. 
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(b). E-plane field patterns for A 2 = 0.75 cm. VSWR = 1.05. 

Figure 32. Moment method and UTD solutions for aperture-matched horn 
as the radius of curvature (A 2 ) of the rolled edge section 
is varied. 
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(c). E-plane field patterns for A 2 = 1.50 cm. VSWR = 1.04. 
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(d). E-plane field patterns for A 2 = 3.00 cm. VSWR = 1.04. 


Figure 32. (Continued). 
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(e) . E-plane field patterns for h? 


7.62 cm. VSWR * 1.04. 
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(f). E-plane field patterns for A 2 = 9.00 cm. y VSWR = 1.04. 


Figure 32. (Continued). 
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(a). Two-dimensional geometry for aperture-matched horn (E-plane), 
f = 10 GHz, w = 1.016 cm, h = 7.85 cm, Aj = 3.05 cm, 

A 2 = 7.62 cm, a varies. 
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(b). E-plane field patterns for a = 15. 


VSWR = 1.07. 


Figure 33. Moment method and UTD solutions for aperture-matched horn 
as the half-flare angle (a) is varied. 
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(d). E-plane field patterns for a = 45. VSWR = 1.04. 
Figure 33. (Continued). 
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(e). E-plane field patterns for a = 60. VSWR = 1.03. 
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(f). E-plane field patterns for a = 75. VSWR = 1.03. 
Figure 33. (Continued). 
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remains relatively constant. The corresponding UTD results are also 
shown in Figure 33. The UTD solutions agree well with the moment method 
patterns for a less that 60°. However, for a = 60°, the UTD results 
begin to disagree with the moment method results. Since the aperture 
width is held constant, the length R|_ of the straight sides of the horn 
decreases as a increases. Also, the curved throat aperture becomes 
closer to the horn aperture. Thus, the infinite wedge approximation 
breaks down again. Note that if both a and h were increased so that Rl 
did not decrease, the UTD solution would remain valid. 

The moment method patterns shown in Figure 34 indicate that as h 
increases, the main beam exhibits a ripple. This is expected since a 
large aperture is likely to have a greater phase variation across it. 
.The VSWR is relatively constant as h varies. The UTD patterns as h is 
varied are also shown in Figure 34. There is excellent agreement with 
the corresponding moment method results, except for h = 3 cm. For such 
small values of h, the length Rl of the straight sides of the horn is 
small so that the infinite wedge approximation is no longer valid. 

As an example, consider the aperture-matched horn geometry shown in 
Figure 35. It represents the E-plane cross-sectional view of a horn 
which was actually designed for use in a compact range. The design 
requirements include a main beam which is relatively flat with less than 
1 dB of variation within a 50° beamwidth. The backlobes must fall to 
approximately 40 dB below the maximum. The frequency range of 
operation is from 8 GHz to 12 GHz. In addition, a low VSWR is 
necessary since the same horn is used to transmit an incident field and 
receive the resulting scattered field. 
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(a). Two-dimensional geometry for aperture-matched horn (E-plane), 
f = 10 GHz, w = 1.016 cm, a = 45°, Aj = 3.05 cm, 

A 2 = 7.62 cm, h varies. 
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(b). E-plane field patterns for h = 3.00 cm. VSWR = 1.04. 

Figure 34. Moment method and UTD solutions for aperture-matched horn as 
the aperture width (h) is varied. 
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(c). E-plane field patterns for h = 7.85 cm. VSWR = 1.04 
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(d). E-plane field patterns for h = 9.00 cm. VSWR = 1.04 


Figure 34. (Continued) 
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(e). E-plane field patterns for h = 12.00 cm. VSWR = 1.04. 
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(f). E-plane field patterns for h 


15.00 cm. VSWR 


1.04. 


Figure 34. (Continued). 
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The calculated and measured patterns are shown in Figure 36. The 
moment method results show excellent agreement with the measured 
patterns. The UTD solutions are within 1 dB of the measurements inside 
the flare angle of the horn, and they are well within 3 dB of the 
measurements out to the broadside direction. The discrepancies in the 
backlobes may be attributed to the inability to accurately model the 
Structures feeding the horn. Also, higher order terms may be needed in 
the diffraction coefficient. In any case, both the UTD and the moment 
method solutions agree well with the measured patterns in the major 
regions of interest. 

The measured VSWR of the horn is shown in Table 1. As indicated 
by the moment method solution, the VSWR is higher at the lower 
frequencies but excellent after about 9.5 GHz. The VSWR had to be 
traded off against the pattern requirements in that the rolled throat 
section was directly affecting the VSWR and main beam pattern. As the 
throat radius increased, the VSWR improved but the main beam pattern had 
too much taper. As a compromise, the throat radius was chosen to 
provide the pattern requirement which was much more critical in the 
design. One should note that a VSWR of 1.2 is typical of most 
traditional horn antennas, so these results still exceed the VSWR 
performance of most horn antennas. 

Clearly, this design does satisfy the necessary specifications. At 
this point, however, an absorbing material was added to the rolled edges 
of the horn to further improve its performance. 


76 



Another property of interest for such a feed antenna is the phase 
of the field. It is very necessary to have a relatively constant 
phase across the main beam of the feed antenna. Figure 37 shows the 
measured phase for the compact range horn with the absorbing material on 
the rolled edges. Note that the phase remains relatively constant as a 
function of angle as well as frequency. This is characteristic of 
aperture-matched antennas since the main beam is dominated by the 
throat term. 

In conclusion, the aperture-matched horn can be designed using 
these analytic techniques to provide uniform patterns across a wide beam 
width with low backlobes, low VSWR, a wide frequency response, and 
constant phase center. 


TABLE 1 


MEASURED VSWR VERSUS FREQUENCY FOR COMPACT RANGE 
APERTURE-MATCHED HORN ANTENNA WITH ABSORBING 
MATERIAL ON ROLLED EDGES 


Frequency (GHz) 

VSWR 

8.0 

1.16 

8.5 

1.12 

9.0 

1.08 

9.5 

1.04 

10.0 

1.02 

10.5 

1.03 

11.0 

1.04 

11.5 

1.03 

12.0 

1.04 
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a. E-plane field patterns for f = 8 GHz. 


Figure 36. Comparison of calculated and measured E-plane field patterns 
for compact range aperture-matched horn. 
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b. E-plane field patterns for f = 9 GHz, 


Figure 36. (Continued). 
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d. E-plane field patterns for f = 11 6Hz. 


Figure 36. (Continued). 
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E-plane field patterns for f = 12 GHz. 


Figure 36. (Continued). 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 

In summary, an aperture-matched horn can be designed to have a 
broader main beam and lower side lobes than standard pyramidal horns. 

In addition, it has a lower VSWR. For certain applications, 
particularly for compact range feed antennas, these characteristics are 
essential. 

Two methods for calculating the E-plane field pattern for such an 
aperture-matched horn were considered. In both approaches, the actual 
three-dimensional horn geometry was approximated by a two-dimensional 
geometry such that any coupling effects between the E- and H-planes are 
neglected. 

Moment method techniques were used to calculate the field by 
solving an integral equation obtained from the reaction concept. The 
VSWR was determined from an internal field. In fact, this method gives 
accurate results provided one uses at least four current samples per 
wavelength. However, at high frequencies, the number of unknowns 
becomes large resulting in intolerably large amounts of computation 
time. 
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The UTD approach presented in this report helped overcome this 
problem. A diffraction coefficient for a source mounted on a pefectly 
conducting planar surface smoothly terminated by a circular cylinder was 
obtained. Using this, coefficient, the field distribution across a 
curved aperture in the throat of the horn was calculated. Next, this 
distribution was expanded in a series of cylindrical modes. The 
resulting series provided the field incident from the throat of the horn 
such that UTD techniques could then be used to calculate the far field 
pattern. This method requires much less computation time than the 
moment method. However, the UTD Is a high frequency technique; it fails 
at low frequencies where the horn dimensions are electrically small. 

At high frequencies, the UTD solution converges asymptotically to the 
exact solution. Specifically, a discrepancy occurs when the radius of 
curvature of the circular cylinder Is less than one half of a 
wavelength. Also, there is a discrepancy when the distance from the 
source to the junction Is less than 1.5 wavelengths. 

Ordinarily, one would expect the UTD solution to be valid at even 
lower frequencies. The disagreement described above suggests that this 
UTD diffraction coefficient may be Improved. Perhaps the addition of 
higher order terms may help the low frequency behavior. In any case, 
the diffraction coefficient obtained here works very well for 
electrically large horns. 

The UTD solution also becomes Inaccurate when the horn geometry 
loses its wedge-shaped appearance. Then the Infinite wedge 
approximation used to obtain the modal expansion of the throat field 
breaks down. This occurs whenever the length Rl of the straight sides 
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of the horn becomes too small. When this happens, the curved throat 
aperture becomes very close to the planar horn aperture; in fact, it may 
even extend beyond the horn aperture. 

Fortunately, the UTD approach works well for the type of geometry 
typically encountered in practice. For low frequencies where the UTD 
fails, the moment method approach is suitable. Together, the two 
methods may be used to efficiently design an aperture-matched horn to 
Operate over a broad range of frequencies. It is suggested that the UTD 
be used to initially design the horn; then, the moment method should be 
used to "fine tune" the result in that it represents a more complete 
solution. On the other hand, one should beware when the moment method 
and UTD do not agree, particularly at high frequencies. It is entirely 
possible that the moment method is being applied in an incomplete sense 
such as a lack of current segments needed to realistically represent the 
horn geometry under test. At low frequencies, the UTD solution may be 
suspect, but then an accurate moment method solution would involve a 
relatively small number of unknowns and therefore little computation 
time. Together, the two methods allow one to quickly and efficiently 
design an aperture-matched horn. 

This analytic approach was used to design an aperture-matched horn 
to feed a compact range measurement system. The requirements for this 
feed antenna were very stringent and necessitated many iterations on the 
computer. However, the desired performance has achieved based on 
experimental verification. In fact, the calculated results agreed very 
well with the prototype measurements. Thus, this design procedure 
has been sucessfully applied in this application. 
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APPENDIX A 


EVALUATION OF THE DIFFRACTION COEFFICIENT FOR A SOURCE 
MOUNTED ON A PERFECTLY CONDUCTING PLANAR SURFACE SMOOTHLY 
TERMINATED BY A CIRCULAR CYLINDER 

Consider only the hard case, i.e., the electric field transverse to 
the z-direction. The diffraction coefficient for a perfectly conducting 
two-dimensional surface with a discontinuity in surface curvature [12] 
is given by 

e -jir/4 C h (5 1 )F(X 1 ) - C h U 2 )F(X 2 ) 

D h “ /2irk cos 4> + cos <(>' (A.l) 
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in which V(t) and W2(t) are the Fock type Airy functions given by 
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Consider the quantity C^(^)F(X^) for the far field (s-*-«) with 
plane wave incidence (s'-*-*) as a^ approaches infinity. Then 

L • (A. 12) 

X -*• «o (A. 13) 

and 

X I - “ . (A. 14) 


For large arguments, the transition function F( X) approaches unity so 
that 


F ( X ) -► 1 (A. 15) 

and 

p ( x l) ♦ 1 . (A. 16) 

For large negative arguments, the Pekeris function has the asymptotic 
form given by 


q (x) 


- 1/2 / 1 x | (1 - ~f) e~^ x Z 12 e J11 ^ 4 
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Thus for the far field with plane wave incidence, one obtains 
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This term may be modified for the general case, i.e., s and s' finite, 

by multiplying by the transition function F(Xj) which takes into account 

the necessary range dependence. Also for grazing incidence ( <t> 1 =0 ) , the 

parameter is reduced by a factor of four to insure continuity at the 

shadow boundary. 

% 

Hence, in the lit region (0 <<j><tt) , the field for a source mounted on 
a perfectly conducting planar surface smoothly terminated by a circular 
cylinder is given by 

e _jks 

H(Q) = zH i (Q J )D h -7=- (A. 19) 

where ^(Qj) is the field incident on the junction at point Qj, s is the 
distance from Qj to the field point, and the appropriate diffraction 
coefficient is 
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(A. 20) 
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The functions F(x) and q*(x) are given in Equations (A. 8) through 
(A. 11). The geometry is shown in Figure 20. 

In the shadow region ( ir<4><2ir) , there is no incident field, and the 
diffracted field is a creeping wave field. It can be described by the 
diffraction coefficient of Pathak, Burnside, and Marhefka [2] but needs 
to be multiplied by an appropriate constant which is chosen to insure 
field continuity at the shadow boundary. Hence, in the shadow region 
( tt< <|><2tt) , the field for a source mounted on a perfectly conducting 


planar surface smoothly terminated by a circular cylinder is given by 


H(Q) = zH 1 (Qj) C T h ~^ 


(A. 28) 


where s^ is the distance from the creeping wave launch point Qj to the 
field point and 
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The geometry is shown in Figure 21. The constant C is the ratio of the 
field H|_( Q) in the lit region to the field H5 ( Q) in the shadow region, 
both evaluated at the shadow boundary SB such that 

ii 
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Clearly, as one approaches the shadow boundary. 
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and hence, the field is continuous. 

The resulting diffraction coefficient for a source mounted on a 
•perfectly conducting planar surface smoothly terminated by a circular 
cylinder is given by 
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where the various parameters are as defined previously. 
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For suitable dimensions, the coefficient DpQ agrees quite well with 
the numerically derived diffraction coefficient of Burnside and Chuang 
[10] as shown in Figures A.l and A. 2. Note that both the magnitude and 
phase of the diffraction coefficients divided by the source distance s' 
are compared as both the cylinder radius a£ and the source distance s' 
are varied. Clearly, the agreement is good for &2 > 0.5X and 
s' > 1.5X. For the dimensions typically encountered in an 
aperture-matched horn design, the diffraction coefficient Dpc is an 
excellent engineering approximation. Horns of smaller dimensions may be 
easily handled with the moment method. 
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(a). Two-dimensional geometry, s' = 3X. 



(b). Magnitude and phase of Dpc/Zs 1 " for a 2 = 0.5X. 
(MM and UTD .) 


Figure A.l. Comparison of UTD diffraction coefficient and numerically 
derived diffraction coefficient of Burnside and Chuang as 
the radius of curvature (a?) of the cylinder is varied. 
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(a). Two-dimensional geometry. a 2 = 3X. 



(b). Magnitude and phase of Ddq//? 1 ’ for s' = U5X. 
(MM and UTD .) 


Figure A. 2. Comparison of UTD diffraction coefficient and numerically 
derived diffraction coefficient of Burnside and Chuang as 
the source distance (s') is varied. 
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(e). 


Magnitude and phase of Dp C //s r for s' = 3. OX. 
(MM and UTD .) 
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